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g ■ Abstract 

In this work, we examine critically the relation between orbital magnetic dipole (scissors mode) 
strength and quadrupole deformation properties. Assuming a simple K = ground state band in 
an even-even nucleus, the quantities Q(2+) (i.e., the static quadrupole moment) and B(E2)q 1 ->2 1 



both are described by a single parameter — the intrinsic quadrupole moment Qo- In the shell model, 
we can operationally define Qo (Static) and Qq{BE2) and see if they are the same. Following a 
brief excursion to the sd shell, we perform calculations in the fp shell. The nuclei we consider 
( ' 6 ' 48 Ti and 48 ' 50 Cr) are far from being perfect rotors, but we find that the calculated ratio 
Qo (Static) /Qo(B E2) is in many cases surprisingly close to one. We also discuss the quadrupole 
collectivity of orbital magnetic dipole transitions. We find that the large orbital B(M1) strength 
in 44 Ti relative to 46 Ti and 8 Ti cannot be explained by simple deformation arguments. 

PACS numbers: 
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I. INTRODUCTION 



In this work we will make a comparison of the shell model and the collective model for 
several quantities that are sensitive to nuclear deformation. These include B(E2) 1 s, static 
quadrupole moments, and orbital magnetic dipole transitions. This will be a theory versus 
theory work. Some experimental results are quoted and serve as anchors for our results, but 
we will not be inhibited by the lack of experimental data in doing these calculations. We 
plan in the near future to make a more extensive theory-experiment comparison. But there 
are holes in the experimental data which must be filled. 

The main thrust of our work will be to understand the relationship of orbital magnetic 
dipole transitions to quadrupole deformations in the nucleus. For example, after the ex- 
perimental discovery in heavy deformed nuclei of the relation between the orbital magnetic 
dipole strength and nuclear deformation Q], there have been many works which relate the 
orbital Ml (scissors mode) strength to electric quadrupole transition rates (-E>(-E2) 0l ^2i)) 
often assuming that they are proportional to each other 0, 0, U, 0, E|. 

First, though, we shall do survey calculations of _B(_E2)'s and static quadrupole moments 
in the fp shell to see how well the shell model relates to the simple rotational model of Bohr 
and Mottelson In this rotational model, the formulae for B(E2)'s and static quadrupole 
moments involve a single parameter — the intrinsic quadrupole moment. These formulae are, 
respectively, 

B(E2) = j^-Ql(B) \(hK20\I 2 K)\ 2 (la) 

107T 

where B and S stand for B(E2) and "Static", respectively. Here Qq is the intrinsic 
quadrupole moment — what we would see in the rotational frame. On the other hand, Q(I) 
is what we measure in the laboratory. In the simple rotational model, Qo(B) is equal to 

Qo(S). 

For the case I± = 0, I 2 = 2, the Clebsch-Gordan coefficient above is 1. For a simple 
K — band in an even-even nucleus, we obtain 

B(E2)^ 2 = ^Ql(B) (2a) 

l07T 
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Q(2 + ) = ~Q (S). (2b) 

Note that the laboratory quadrupole moment has the opposite sign of the intrinsic 
quadrupole moment — a well known result. It can be understood physically by imagin- 
ing rotating a cigar (which has a positive quadrupole moment, i.e., prolate) about an axis 
perpendicular to the line of the cigar. This will trace out a flat pancake shape which is 
oblate. 

We then find that the ratio 



%™ = - V^^EL = -1.10387054EL (3) 
Q (B) 2 V 16tt ^/bJeT) ^/bJE2) V ' 

II. QUADRUPOLE PROPERTIES IN THE sd SHELL 

Although we will be performing calculations in the fp shell, we shall here briefly look 
over the experimental situation in the sd shell. In Table |J we show experimental values of 
Q(2+) Q, B(E2) Q, and the ratio \Q (S)/Q (B)\ as given by Eq. ©. We also show the 
experimental values of E(4±) / E(2^) as a measure of how close we are to the rotational limit 
of 10/3, or the vibrational limit of 2/1. 

TABLE I: Experimental data on Q(2f) and B(E2) in the sd shell and the ratio \Q (S)/Q (B)\; 
we also give the experimental values of E(A) / E{2). 





0(2f) 
[e fm 2 ] 


B(E2) 
[e 2 fm 4 ] 


\Qo(S)/Q (B)\ 


E(A+)/E(2f) 


20 Ne 


-23 


340 


1.377 


2.600 


22 Ne 


-19 


230 


1.383 


2.634 


24 Mg 


-16.6 


432 


0.802 


3.012 


28 Si 


+16.5 


320 


1.018 


2.596 


32 S 


-14.9 


300 


0.950 


2.000 


36 Ar 


+11 


340 


0.658 


2.240 


40 Ar* 


+ 1 


330 


0.061 


1.980 
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Note that the static quadrupole moments of 20 Ne, 22 Ne, 24 Mg, and 32 S are negative, while 
those of 28 Si and 36 Ar are positive. If we limit ourselves to axial symmetry, this indicates 
that the first group has prolate ground state bands and the second group has oblate ones. 
Skyrme II Hartree-Fock results by Jaqaman and Zamick |10J] correctly give the signs of all 
the static quadrupole moments. The small static quadrupole moment of 40 Ar is consistent 
with magnetic moment results of the 2f state by Stefanova et al. 

The ratio \Qq(S)/Qq(B)\ for 20 Ne is larger than the rotational limit, 1.377 versus 1; 
likewise 22 Ne. In the case of 24 Mg, \Q (S)/Q (B)\ is smaller than for 20 Ne or 22 Ne, despite 
the fact that the spectrum is closer to rotational for 24 Mg. Also surprisingly for 32 S, the 
ratio E(A)/E(2) is 2.000, the vibrational limit, for which one might expect a near zero static 
quadrupole moment. But the ratio \Q (S) / Q (B)\ is 0.950, close to the simple rotational 
prediction of unity. 

In general, it is difficult to correlate \Q (S) / Q (B)\ with E(A)/E{2) assuming a simple 
axially symmetric rotor. 

We should mention that an analysis of the relationship of Qo{S) and Qn(B) has already 
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and Bender et al. I13L albeit not for 



been performed by Bender, Flocard, and Heenen 
the /p-shell nuclei considered here and using a different method. They perform angular 
momentum projections on BCS-Hartree-Fock states obtained with the Skyrme interaction 
SLy6 for the particle-hole channel and a density-dependent contact force in the pairing 
Their calculations are mainly in the sd shell Q| and neutron-deficient lead 
13| . For one nucleus in common, 40 Ca, their results for Op — 0/i, 2p — 2h, Ap — Ah, 
Qp — 6h, 8p — 8h, and 12p — 12h do not differ so much from previous calculations of Zheng, 
Berdichevsky, and Zamick |l4| as far as the intrinsic properties are concerned, but their 
calculation has the added feature of providing an energy spectrum and expectation values 
in the laboratory frame. 

In Ref. the authors predict that 36 Ar is oblate. This is confirmed by the fact that 



channe 
region 



the static quadrupole moment of the 2t state is positive: +11 e fm 2 The experimental 
B(E2) is 340 e 2 fm 4 and |/3 2 | = 0.273 [9]. Using Eq. ©, we find 

Qo(S) 



Qo(B) 

The energy ratio is 



0.6505236 . (4) 



E(At) 4414.36 

— r-TT = = 2.240 . 5 

E{2+) 1970.35 V ; 



These results are consistent with a nucleus not being too rotational. 
The responding nu.be, in the cation of Bender et a.. Q are 

Q(2+) lab = 13 e fm 2 , B(E2) |= 220 e 2 fm 4 , (3 = -0.21 . (6) 

The calculated ratios are 

Q °( S ) - n 9675 E ^ - 2 6545 (7) 

Q (B) ~ ° J675 ' 25(2+) ~ 2 - 6M5 ' {7) 

These calculations (l2| give a more rotational picture than experiment. There is a consis- 
tency, however, in that a larger ratio E{A)/E(2) yields a larger ratio Q (S)/Q (B). 

III. SHELL MODEL CALCULATIONS OF B(E2) AND Q(2+), AND HOW THEY 
RELATE TO THE SIMPLE ROTATIONAL MODEL 

We will put the above relation (jSJ) to the test in a shell model approach for the following 
nuclei: 44 Ti^Ti, 48 Ti, 48 Cr, and 50 Cr. We use the OXBASH program [15] and the FPD6 
interaction 1 6 . 

The nuclei that we have chosen are far from being perfect rotors. Their description falls 
somewhere between vibrational and rotational. The ratios E(4)/ E(2), which would all be 
10/3 in the simple rotational case, are as follows: 1.922, 2.010, 2.118, 2.459, 2.342, for 44 Ti, 
46 Ti, 48 Ti, 48 Cr, and 50 Cr, respectively. 

We perform shell model calculations in a complete fp space using the FPD6 interaction. 
We assign effective charges of 1.5 for the protons and 0.5 for the neutrons. We calculate 
B(E2) 0l ^2 1 and Q{2 + ) (the laboratory Q, of course) and put them into Eq. (jHJ) in order to 
get operational values of Qo(S)/Qo(B). The results are given in Table ITT1 

Except for 48 Ti, the FPD6 results for the ratios are all greater than 0.9, reading a maxi- 
mum of 0.9892 for 48 Cr. It is somewhat surprising that these ratios are so close to 1, given 
that the ratios E(4)/E(2) are much further away from the rotational limit 10/3. 

We can also obtain some of the above ratios from experiment. We refer to the compilation 
of nuclear moments of Stone [8] and of B(E2Ys by Raman et al. j^. Taking these experiments 
at face value, we see that the ratio Qo(S)/Qo(B) reduces to about 0.75 for 46 Ti and 48 Ti, 
but is bigger than 1 for 50 Cr. 

It should be noted that in the simplest version of the vibrational mode, Qo(S) is zero. We 
can imagine a nucleus vibrating between a prolate shape and an oblate shape, and causing 



TABLE II: The quantity Qq{S) / Qq{B) as obtained in the shell model with the FPD6 interaction 
and experiment. 





0(2f) 

V I ) 

[e fm 2 ] 


B(E2) 
[e 2 fm 4 ] 


On(S)/Oa(B) 


Theory (FPD6V 

44 Ti 


-20.156 


607.24 


0.9029 


46 Ti 


—22.071 


682 06 

UU^j • \J\J 


9329 


48 Ti 


-17.714 


560.78 


0.8257 


48 Cr 


-33.271 


1378.4 


9892 


50 Cr 


-30.955 


1219.0 


0.9787 


Experiment: 

44 Xi 




650 




46 Xi 


-21 


950 


0.7521 


48 Ti 


-17.7 


720 


0.7282 


48 Cr 




1360 




50 Cr 


-36 


1080 


1.2092 



the quadrupole moment to average to zero. On the other hand, the B(E2) 0l ^ 2l is quite large 
in this vibrational limit, causing the ratio Qo(S) / Qo(B) to be zero or, in more sophisticated 
vibrational models, quite small. 

IV. RESULTS WITH AN ALTERNATE INTERACTION T0FPD6 

For systems of identical particles, e.g., the tin isotopes, which in the simple shell model 
involve only valence neutrons, one does not get rotational behaviour. One does go closer 
to the rotational limit when one has many open-shell neutrons and protons. Whereas two 
identical nucleons must have isospin 1, a neutron and a proton can have both isospin and 
1. This suggests that the T = part of the nucleon-nucleon interaction plays an important 
role in enhancing nuclear rotational collectivity. 

In this section, we will use an interaction TOFPD6 that is the same as the FPD6 inter- 
action for T = 1 states, but vanishes for T = states. We thus expect that the rotational 
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collectivity will be reduced. This interaction serves as a counterpoint of the full interaction 
in the previous section. It has been discussed before by Robinson and Zamick 2l| . 
We present the results for T0FPD6 in Table ITTT1 

TABLE III: The same as Table [HI but with the T = two-body matrix elements of FPD6 set to 
zero (T0FPD6). 



Nucleus 


QVt) 

[e fm 2 ] 


B(E2) 
[e 2 fm 4 ] 


Qo(S)/Q (B) 


44 Xi 


-0.880 


375.09 


0.0502 


46 Ti 


-8.195 


432.81 


0.4348 


48 Ti 


-8.777 


401.97 


0.4832 


48 Cr 


-21.437 


813.06 


0.8299 


50 Cr 


-20.985 


736.60 


0.8535 



We see that both Q(2f) and B(E2) decrease in magnitude, consistent with the above 
discussion. However, Qo(S) decreases more rapidly than ^JB(E2), so the ratio Q (S)/Q (B) 
is less for this case than when the full interaction is present. For 44 Ti this ratio is 0.0502, 
close to the vibrational limit of zero. For 48 Cr the ratio decreases from 0.9892 to 0.8299; 
nevertheless, it is still substantial, indicating that the T = 1 interaction, acting alone, can 
lead us to some extent in the direction of the rotational limit. 



V. RANDOM INTERACTION STUDIES 



Nuclei in the region we are conside ring have undergone Random Interaction studies. We 



refer to the works of Velazquez et al. 17( and Zelevinsky and Volya [18[ . These works were 
stimulated by that of Johnson, Bertsch, and Dean 

In particular, in the work of Zelevinsky, a quantity is considered which is proportional 
to the square of Qq{S)/Qq{B) and which is normalized to (2/7) 2 if Q (S)/Q (B) equals 1. 
He calls this the Alaga ratio. He selects cases for which the random interaction yields a 
J = 0, J = 2 sequence of lowest energies. For these he finds two peaks, one corresponding 
to Qo(S)/Qo(B) = 1 and the other to Q (S)/Q (B) = 0. In our terminology, these would 
correspond to the rotational limit in the former case and either the simple vibrational limit 
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or the spherical limit in the latter case. He cites early work of R. Rockmore [20] as affording 
an explanation of this surprising behaviour. 

VI. ORBITAL MAGNETIC DIPOLE TRANSITIONS IN 44 Ti, 46 Ti, AND 48 Ti 
The orbital magnetic isovector dipole transitions, i.e., scissors mode excitations, also 

n 

display collective behaviour There are systematics which suggest that i?(Ml) orbital is 
roughly proportional to B(E2). There are more detailed sophisticated relations as well. If 
one uses a simple quadrupole-quadrupole interaction, the energy weighted 5(Ml) orbital is 
proportional to the difference ^i?(E2)i SOSCa i ar - B(E2) iBovectoT J 3J]. 
The bare orbital Ml operator is 

y^EWaCO, (8) 

where gi is 1 for a proton and for a neutron. This is the operator that we use in the 
calculations. 

How to extract the scissors mode strength is not completely unambiguous. The mode is 
associated with low-lying 1 + excitations at around 3 MeV. But the strength can be frag- 
mented even at this lowest energy. Besides this, there is orbital strength at higher energies, 
a somewhat grassy behaviour where individual states are very weakly excited but, because 
there are so many of them, the total orbital strength can be significant. 

Therefore, we will give three sets of values (see Table ITVj) . First, we give the strength to 
the lowest state, then to the lowest 10 states, and finally to the lowest 1000 states (except 
for 48 Ca, where we include only 300 states). The 10-states strength should encompass what 
we usually call the scissors mode, while the 1000-states strength is close to the total strength 
including the grassy, non-collective part. It would appear that the highest excitation energies 



reached in the experiments ^,|22j are not sufficient to reach the T + 1 part of the spectrum. 

We first discuss the nuclei 46 Ti and 48 Ti, for which there are some data on B(M1). We see 
consistently that the orbital B(M1) strength is larger in 46 Ti than in 48 Ti. This is consistent 
with the fact that 46 Ti has a greater B(E2) and static 2 + quadrupole moment than 48 Ti. 

We next consider the N = Z nucleus 44 Ti, for which there is no data because this nucleus 
is unstable. The isoscalar orbital B(M1) strength is very weak. This is also true for the 
spin B(M1), but for a different reason. The isoscalar spin coupling is much smaller than 
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TABLE IV: The calculated orbital Ml strengths (fJ^)- Unless indicated, the calculations are made 
with the FPD6 interaction. 



44rp; 46r-n; 48^; 48rn; /46rp; 48 



T -> T 



Lowest state 


0.0017 


0.305 


0.105 


0.3443 


Lowest 10 states 


0.0320 


0.5979 


0.3056 


0.5111 


Lowest 100 states 




0.79 


0.504 


0.6380 


All states 


0.0355 


0.9195 


0.7191 


0.7820 


All states (T0FPD6) 


0.0583 


0.2166 


0.2235 


1.0319 



T -> T + 1 

Lowest state 0.862 0.0991 0.0041 0.4450 0.784 

Lowest 10 states 1.4317 0.368 0.1951 0.5302 1.3855 

Lowest 100 states 2.12 1.994 

All states 2.127 0.5616 0.3099 0.5518 2.271 a 

All states (T0FPD6) 1.1118 0.3940 0.1671 0.4241 

"Lowest 300 states 

the isovector one. For the orbital case, the couplings are equal because the operator is 
S protons ^- F° r the orbital case, the B(M1) isoscalar is very weak because the correlations 
due to the nuclear interaction move the ground state towards the 577(4) limit, in which LS 
coupling holds and for which the ground state is a pure L = state. For this extreme case, 
the B(M1) orbital isoscalar will vanish. 

The transitions of interest for 44 Ti are, therefore, the isovector orbital dipole ones. The 
(T — ► T + 1) B{ Ml) orbital summed strength is larger in 44 Ti than the (T — > T) and [(T — > 
T) + (T -»• T+l)] strengths in 46 Ti and 48 Ti, which are 0.5615 (1.4810) and 0.3099 (1.0288) 
fiff, respectively. On the other hand, 44 Ti is not more deformed than 46 Ti. According to 
Raman et al. [9J, the values of the quadrupole deformation parameters (3 for 44 > 46 > 48 Ti and 
48,50Q r are ^ reS p ec tively, 0.27, 0.317, 0.269, 0.335, and 0.293. Thus, we have here in fp- 
shell nuclei a counter-example to the experimentally established proportionality between 
the orbital B(M1) and the B(E2) in heavy deformed nuclei. Perhaps there are correlations 
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which cause an enhancement for N = 2 nuclei. 

An analysis by Retamosa et al. [23[ in the sd shell comparing 20 Ne, 22 Ne, and 24 Mg was 
performed (somewhat analogous to 44 Ti and 46 Ti for the first two cases), but no anomaly was 
reported there. In the SU(3) model, they found consistency in the relation of -B(Ml) or b>itai 
to deformation. In this limit, the summed Ml strengths (all orbital) for 20 Ne 22 Ne, and 
24 Mg were 1.1, 1.17, and 1.6 /z^, respectively. Looking at the Raman tables |9| for these 
nuclei, there is some complication — the deformation parameters fa are not in one-to-one 
correspondence with the B(E2)'s. The values of (B(E2), fa for these three nuclei from 
the Raman tables jj are, respectively, (0.034,0.728), (0.0236,0.562), and (0.0432,0.606), 
where the units for B(E2) are b 2 . The authors also do calculations with a more realistic 
interaction, but no enough strengths are listed in order to make a comparison for the point 
we are trying to make. Retamosa et al. [23] also give strengths to the first 10 + states in 
44 Ti; our numbers are consistent with theirs. 

Earlier works on the shell model for light nuclei include L. Zamick j^j and A. Poves 



VII. CLOSING REMARKS 



In this work we have examined what predictions the shell model make for collective 
properties which are after dealt with in the rotational model. Although the nuclei are far 
from perfect rotors, the calculated ratio Qo(S)/Qq(B) is fairly close to 1 in many cases. 
When the FPD6 interaction is used, the orbital magnetic dipole transitions for 46 ' 48 Ti also 
fit into this picture, although there is the added complication of separating the collective 
from the non- collective part in this case. Also there is a substantial enhancement for the 
N = Z nucleus 44 Ti, which cannot be explained as purely a deformation effect. We hope our 
work will stimulate more experimental investigations. There is information of B (Ml) rates 
in 46 Ti and 48 Ti, but thus far the orbital B(M1) has only been extracted in 48 Ti. However, 
in a short time, we will be able to make a more extensive theory-experiment study of these 
magnetic dipole transitions. 
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